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Domain-wall solutions in four-dimensional supersymmetric field theories with 
distinct discrete vacuum states lead to the spontaneous breaking of super- 
symmetry, either completely or partially. We consider in detail the case 
when the domain walls are the BPS-saturated states, and 1/2 of supersym¬ 
metry is preserved. Several useful criteria that relate the preservation of 
1/2 of supersymmetry on the domain walls to the central extension appear¬ 
ing in the N = 1 superalgebras are established. We explain how the cen¬ 
tral extension can appear in V = 1 supersymmetry and explicitly obtain 
the central charge in various models: the generalized Wess-Zumino models, 
and supersymmetric Yang-Mills theories with or without matter. The BPS- 
saturated domain walls satisfy the first-order differential equations which we 
call the creek equations, since they formally coincide with the (complexi¬ 
fied) equations of motion of an analog high-viscosity fluid on a profile which 
is given by the superpotential of the original problem. Some possible ap¬ 
plications are considered. We also briefly discuss BPS-saturated strings. 



I. INTRODUCTION 


Recently a new mechanism for supersymmetry (SUSY) breaking was discussed in 
The essence of the idea is that higher dimensional theory is compactihed on a topological 
(soliton-like) defect, not invariant under supersymmetry transformations. In the theories 
where several discrete vacua arise, the domain wall solutions take place. Since each given 
solution of this type breaks a part of the translational invariance of the theory, it is quite 
natural that in the sector with the domain wall supersymmetry is spontaneously broken too. 
In Refs. it was observed that in some instances a part of supersymmetry is preserved. 
The simplest example is the Wess-Zumino model [^. The domain wall in this model was 
shown [0] to be the Bogomolny-Prasad-Sommerfeld (BPS) 0 saturated state, thus general¬ 
izing the old result referring to kinks in two-dimensional theories. In the present paper 
we analyze the BPS-saturated domain walls in various models in more detail. The very 
existence of such states is associated with a central extension of Y = 1 superalgebra. Since 
the possibility of such central extensions seemingly contradicts the well-known theorems, we 
hrst explain (Sect. 2) what type of central extension can arise, and how the general asser¬ 
tions are circumvented. The key point is that the central extension is automatically zero for 
all spatially localized held conhgurations. It need not necessarily vanish, however, for those 
held conhgurations that interpolate between distinct vacua at spatial inhnities, i.e. in the 
presence of the domain walls. The concrete form of the centrally extended N = 1 superal¬ 
gebra is derived for the generalized Wess-Zumino models. Here the central charge appears 
at the classical (tree) level. The central charge in the gauge theories without superpotential 
emerges at the one-loop level, as a quantum anomaly. We give an explicit diagrammatic 
derivation of this anomaly for supersymmetric QED (SQED). We then comment on how 
this anomaly is related to those previously established in the non-Abelian theories |p-|^ 
and present the centrally extended N = 1 superalgebra for generic non-Abelian theory with 
matter and superpotential. 

In Sect. 3 we return to the generalized Wess-Zumino model, and derive the classical equa- 
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tions defining the BPS-saturated domain walls. These equations turned out to be analogous 
to those describing the mechanical motion of a high viscosity fluid on a multidimensional 
prohle with mountain ridges and valleys (in the minimal Wess-Zumino model the prohle is 
two-dimensional). A rich physical intuition one has in mechanical motion helps us establish 
conditions for the existence of the solution(s). If the theory has more than one discrete 
vacuum state the domain wall as the solution of the classical equations of motion always 
exists. The BPS saturation equations are stronger constraints, which may or may not have 
solutions. We establish several useful criteria for the solutions to exist. 

The last part of this section is devoted to BPS-saturated strings. The central charge 
need not vanish in the case of string-like topological defects. If it is non-zero, there arises a 
possibility of having stable string-like solutions preserving 1/2 of SUSY. We address some 
general aspects of the BPS-saturated strings and then consider a specihc example, in a 
generalized Wess-Zumino model. 

Section 4 is devoted to the minimal Wess-Zumino model. After the domain wall is formed 
one can view the original (3 -|- 1) dimensional theory as a low-energy (2 -|- 1) dimensional 
reduction. We build the reduction explicitly. All helds of the original theory are decomposed 
in the modes in the wall background. A special role belongs to zero modes - the only ones 
surviving in the limit of vanishing energy. These modes are localized on the wall. Among 
the non-zero modes some are localized, others are not. We construct (2 -|- 1) dimensional 
superhelds which realize the residual N = 1 SUSY in three dimensions (this is Y = | SUSY, 
from the point of view of the original four-dimensional theory). We obtain a Lagrangian with 
an inhnite number of interacting superhelds, in (2 -|- 1) dimensions. The resulting theory is 
explicitly supersymmetric and thus have vanishing corrections to the vacuum energy density. 
In terms of the four-dimensional theory this means that the domain wall energy density is 
not renormalized. The non-renormalization theorem for the energy density e in the BPS- 
saturated domain walls was established in []^,^. This phenomenon is related to the existence 
of the central extension of Y = 1 superalgebra, much in the same way as the exact expression 
for the monopole mass in the Y = 2 SUSY Yang- Mills theory can be traced back to the 
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central charge appearing in this theory |P . 

The representation of (2 + l)-dimensional supersymmetry is minimal, since there is no 
analog of the chiral supertransformation in (2 + 1) dimensions. Unlike the situation in 
four dimensions, superspace can not be split into two chiral subspaces. This circumstance 
explains why there are no non-renormalization theorems for superpotential of the type we 
are used to in four dimensions |^. As a result, in the (2 -|- 1) dimensional theory with = 1 
SUSY, the shape of the domain wall is subject to perturbative corrections. 

In Sect. 4 we also calculate the one-loop perturbative correction to the prohle of the 
domain wall in the minimal Wess-Zumino model. 

Section 5 is devoted to sample applications of the ideas developed in this work. First 
we consider the tunneling problem in multidimensional (non-supersymmetric) quantum me¬ 
chanics with two or more degenerate classical minima. For a class of potentials of a special 
form the tunneling trajectory (in the Euclidean time) is determined by the same hrst-order 
equations as the domain wall prohles in the generalized Wess-Zumino models. The fact 
that these equations do have solutions can be established from the global features of the 
potentials, without actual solutions. The analytic expression for the classical action on these 
least-action trajectories follows immediately. This classical action determines the tunneling 
exponent. 

In the second part of Sect. 5 we make use of the calculated quantum corrections to the 
prohle of the domain wall. They are related to the subleading corrections in the asymp¬ 
totic behavior of the multiparticle production amplitudes at threshold when the number of 
particles tends to inhnity |1l0|,pd|. By exploiting this relation we hnd the hrst subleading 
correction to the corresponding asymptotic formula in the minimal Wess-Zumino model. 


This correction was previously known in non-supersymmetric theories 11,13,12 


In Sect. 6 we address the issue of renormalization of the soliton mass /i in the Witten- 
Olive solution of the (1 -|- 1) dimensional theory. The reason why /i is renormalized while £ 
in four-dimensional theories is not, is revealed. 

Technical details of the formalism extensively exploited throughout the paper are col- 
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lected in Appendices. 


Some of the results referring to the generalized Wess-Zumino model, to be considered 


below, were discussed previously, mainly in the stringy context, see and references 

therein. In particular, the arguments at the end of Sect. 3.A can be found in these works. 
We present them for the sake of completeness. Our main emphasis is put on new results: 
(i) the anomalous term due to the gauge fields appearing in the central extension of A^ = 1 
superalgebra; (ii) criteria for the existence of the BPS-saturated domain walls which is far 
from being guaranteed by the occurrence of a non-vanishing central extension; (iii) so far 
unexplored applications. 


II. CENTRAL EXTENSIONS 


A. Generalities 


In this section we explain how and what kind of central extension can appear in A = 1 
super algebra 


{QqQ/b} — ‘2Ta/3 , 


( 1 ) 


where Qa are the supercharges, a = 1,2 (the Hermitian conjugated supercharges will be 
denoted Qa)- It is common wisdom, that N = 1 supersymmetric theories can not contain 
central extensions. The general classification of superalgebras dates back to the classical 
paper [|T^]. Should the central extension appear in the anticommutator ([^), it will clearly 
belong to (0,1) representation of the Lorentz group. Here Qa is the supercharge, and 
spinorial notation is used throughout the paper, if not stated to the contrary. This fact 
is obvious, since the central extension should be symmetric with respect to two undotted 
indices. The existence of an extra conserved quantity that is not a Lorentz scalar, in addition 


to four-momentum, is forbidden by the Coleman-Mandula theorem for all Lorentz- 
invariant theories with non-trivial S matrix. Thus, in order to have a central extension we 
must violate some assumptions of this theorem. Let us look at the theories where the Lorentz 
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symmetry is spontaneously broken. This situation can be realized in the theories with 
spontaneously broken discrete symmetries, with two or more discrete degenerate vacua. Field 
configurations that interpolate between the distinct vacua are extended objects which are 
not invariant under the action of the Poincare group. If we choose such a field configuration 
as our “vacuum”, then we are clearly in the situation with the spontaneously broken Lorentz 
symmetry. In this case the existence of a central extension does not contradict the Coleman- 
Mandula theorem, Tal3 7^ 0. The central charge, however, must vanish in the sector with 
the Lorentz-invariant vacuum. 

In supersymmetric theories one can have a discrete set of (degenerate) vacua even without 
spontaneous breaking of any discrete symmetry. Supersymmetry, rather than the sponta¬ 
neous breaking of the discrete symmetry, will keep the energy of these states at zero, so that 
the interpolating field configurations, the domain walls, can be stable both classically and 
quantum-mechanically. As a matter of fact, this situation is quite typical. Supersymmetric 
gluodynamics, the simplest SUSY gauge theory, does have N + 1 vacuum states (for the 
gauge group SU{N)) that are not continuously connected. N vacua are chirally asymmet¬ 
ric and are related to the spontaneous breaking of the Z 2 N symmetry; one extra vacuum 
is a recently discovered chirally symmetric state whose existence is unrelated to the 
spontaneous breaking of any discrete symmetry. In the generalized Wess-Zumino models 
multiple vacuum states unrelated to the spontaneous breaking of any discrete symmetry are 
quite conventional too. The BPS-saturated domain walls are abundant in these theories, 
the phenomenon is not at all exotic. 

The supercharges appearing in the superalgebra (|l|) are defined as 


Qa = j d^xJl- 


( 2 ) 


where is the conserved supercurrent of the supersymmetric model at hand, in the spin- 
vector form. 




( 3 ) 
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where is the supercurrent in the spinorial form, and the matrices a are dehned in 
Appendix A, where all our notations and conventions are collected. The current conservation 
can be written as 



( 4 ) 


The supercurrent is a member of the supercurrent supermultiplet 



( 5 ) 


where is the Rq current and is the energy-momentum tensor . The relation 


between and is given in Appendix A, Eg. ( |A41D . 

With these dehnitions in hands we can proceed now to a detailed discussion of the issue 
how Tap 7 ^ 0 appears. 


B. Wess-Zumino Model 


In the non-gauge theories, the presence of the central extension can be seen at the tree 
level. The most clear-cut example is the Wess-Zumino model, where we can explicitly End 
the solution of the classical equations of motion interpolating between two vacua of the 
model (the domain wall). In this case the central extension appears already at the tree 
level. The Wess-Zumino Lagrangian in terms of superhelds has the following form. 



( 6 ) 



( 7 ) 


where the parameters p and A can be always chosen to be real. If p = 0, the Lagrangian 


is invariant under the i?o"rotations |^ , 



( 8 ) 
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at the classical level. If /x 7 ^ 0 this invariance is gone, but a discrete subgroup persists, 


0 , ^ i'^ 


F. 


(9) 


This discrete part of Rq symmetry does not mix the real and imaginary parts of the bosonic 
held. The Z 2 symmetry is spontaneously broken, and the corresponding domain wall solution 
(see Sect. 4), that interpolates between two degenerate vacua, exists. The supercurrent has 
the following form (we write everything in the chiral notation, see Appendix A), 


= 2V2 


1 r 




( 10 ) 


Its dependence on the superpotential comes through the equation of motion for the F term. 


F 


dW 
d(j)+ ' 


( 11 ) 


The second square brackets in Eq. ( 0 ) contain a full spatial derivative and, in principle, 
could have been omitted. Upon inspecting Eq. (p! 0 D it becomes clear that the anticommutator 
of two supercharges does not vanish, and is completely determined by the canonic equal-time 
anticommutator 


{V’a(a:)V’a(?/)}e.t. = laaS{x-y). 


( 12 ) 


The term with 9 o 0 ’'"(c’'o)/ 3 o drops out since it can not yield symmetric in (a,/9) structures, 
and then we get for the anticommutator of two supercharges. 


{Qa Qf}} = (-4f)(<? 


RxV m - -<!> 




1 


(13) 


The matrix aap is dehned in Eq. ( |A12|) and is automatically symmetric in a, (3. It converts 
the spinorial indices of the representation (1,0) of the Lorentz group in the vectorial form. 
The expression on the right-hand side is formally ambiguous. Indeed, let us drop the second 
term in the supercurrent ([T0|). This term is a full spatial derivative and should have no 
impact on the supercharge. However, if we calculated {QaQ/s}, using only the hrst term in 
Eq.(pTD for dehning the supercharge Qa, we would get 
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{Qa Q/3} = (-4*)(a)„/3 / d^xVW |,-=o, 


(14) 


instead of Eg. (|T^) . This ambiguity does not affect the value of the central charge, however. 
The term $c}VV/55> is a full superderivative, and it vanishes in any supersymmetric vacuum. 
Thus, whatever wall is considered, /d^ gives no contribution in Eg. (p!3[) . 

We pause here to make a remark regarding the impact of loop corrections in Eg. (p!3D. 
Strictly speaking, when one includes the guantum loops, the tree-level result (P!B|) acguires 
an additional term on the right-hand side, due to the guantum anomaly. This term is 

-7D^(|)<h) 


where 7 is the anomalous dimension of the superheld <h, cf. the last term in Eg. (^ 0 ]) . As 
we already know, such terms (full superderivatives) can be freely omitted. 

Summarizing, the expression for the central charge given in Eg. (|T^) is exact. The operator 
on the right-hand side, being sandwiched between any states in the sector with the unbroken 
Lorentz symmetry (no domain walls or other topological defects) vanishes, in accordance 
with the Coleman-Mandula theorem. In the sector with the domain wall, however, the 
expression is non-zero. It is saturated by the extended held conhguration of the domain wall 
and is proportional to the difference between the vacuum expectation values of W in two 
distinct vacua, between which the given domain wall interpolates. 

In the example considered the value of TV in both vacuum states is real. This is not 
necessarily the generic case. If one limits oneself to consideration of the given pair of vacua 
(or the model considered has exactly two vacuum states) it is always possible to ensure that 
W in both vacuum states is real. Indeed, hrst we can use the fact that the superpotential 
is dehned up to an additive constant. Using this freedom one can always choose >V*i = 0. 
Here the subscript *1 marks the value of the superpotential in the hrst vacuum. Moreover, 
the superpotential is dehned up to a phase. The phase of W can be always rotated by an 
appropriate rotation of the Grassmann variable 6. Using this freedom it is always possible 
to make kV *2 real. It is even possible to make kV *2 real and positive. Below this choice 
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will be referred to as standard, and in many instances we will exploit it. (Of course, if one 
discusses simultaneously two or more distinct domain walls in the problems with three or 
more vacua it may not be always possible to kill the phase.) 

It is convenient to define the central charge under discussion in the following general 
form. Let the charges be normalized in the standard way, 

{QaQp} = ‘^Pa/3y 

and assume that the wall lies in the xy plane (this can always be achieved by an appropriate 
choice of the reference frame). If 

{QaQp] = (-2i)(ri)„/3SA, (16) 

the parameter S will be referred to as the central charge. Here A is the area of the wall. 
From the consideration above it is clear that in the Wess-Zumino model 


S 


2 


W{z = +cxo) — W{z 



(17) 


For the superpotential ( 0 ) the vacua are at 

0*1,2 = ± 

and 

Although the derivation was given for the minimal Wess-Zumino model, Eqs. (|^ and 
(O) are obviously valid for any number of the chiral superfields and arbitrary superpotential. 

Since the supercharges are conserved, the right-hand side of Eq. (^) is a conserved quan¬ 
tity too. This observation can be viewed as an alternative proof of the absence of renormal¬ 
ization of the superpotential. It is well-known [^j that superpotentials are not renormalized 
perturbatively. Non-perturbative renormalization can occur only if Eq.(^) has a non-trivial 
anomaly (see below). 
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The commutator relation (pTl) contains only Q’s not Q’s. This means that the central 
charge S must be an analytical function of various parameters characterizing the problem 
(masses, coupling constants, and so on). 

The fact that we were able to hnd an expression for the central extension as a local 
operator does not necessarily mean that the central charge is non-vanishing. For instance, 
if the parameter /i in Eq.(]^) is put to zero, E = 0. In Ref. |l| examples are given of the 
generalized Wess-Zumino models with two classically distinct vacua such that 

>V(z = -l-cxo) = W^z = —cxo). (18) 


In these models S vanishes too. The domain walls are still present, of course, but they are 
not BPS-saturated. 


Assume now that S 7 ^ 0. By combining Eqs. (|T^) and ([T 6 |) we can get a bound for the 
domain wall energy density e. 

Theorem 1. e > |E|. If the inequality is saturated, i.e. e = |E|, 1/2 of SUSY is 
unbroken on the given domain wall. 

To prove the theorem consider a linear combination of operators Q and Q, 


= (19) 

where (3 is an arbitrary complex parameter, with two components, (3 = {/d^,/3^}. Both 
are treated as c-numbers rather than the Grassmann numbers. Since K is Hermitian, the 
expectation value over any state S 

{S\K^\S) > 0. (20) 


Equation (pO]) implies, in turn, that for any f3 

S{\f3^\^ + l/ST) + 2 Im A > 0 ( 21 ) 

where £ is the energy of the state S. To get the best possible bound on E we must optimize 
our choice of (3. If 
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S = 


(22) 


optimization is evidently achieved at (3^ = —ij3'^ = where p and b are positive numbers, 
while a is a phase. Equation (^2]) then reduces to 

E>pA (23) 


or, in other words. 


e>|S|, (24) 

q.e.d. The walls for which this inequality is saturated (i.e £ = |S|) are called BPS-saturated. 
It is clear that the BPS saturation can be achieved only provided that a linear combina¬ 
tion of the supercharges, acting on the wall, annihilates it. This means, that a part of 
supersymmetry is preserved in the sector with the given domain wall. 


C. Central charge in = 1 gauge theories. 


Let us discuss now supersymmetric gauge theories with or without the superpotential. As 
a warm-up exercise we consider supersymmetric electrodynamics (SQED). The Lagrangian 

is 

C = j (feW^ + H.c)j + \J {Se^S + f e’^T) + j d^OS T + H.c)j (25) 

where S and T are chiral superfields of the electric charge -|-1 and —1, respectively. The 
Lagrangian (^5] ) describes photon, electron and two selectrons. The expression for the 
supercharge can be read off Eq. (|A45|) . 

The model possesses the Rq symmetry 


A ^ e*"A , e , 0, 


S.t 




(26) 


which is broken down to Z 2 by the mass term and the quantum anomaly. The vacuum is 
non-degenerate and, thus, we do not expect to find the domain walls. Nevertheless, it is still 
useful to complete the calculation of Tap. The supercurrent of the model is 
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(27) 


'^a/3/9 — 2 2i Gf}QiX^ QSjjaDX^ + a/ 2 i 


Terms with the total spatial derivatives in the supercurrent, see Eq.( A45 ), are unimportant 
for the present discussion and we omitted them from the beginning. In the tree approxima¬ 
tion the anticommutator {QaQis} is the same as in the Wess-Zumino model with 


W = moS T, 


(28) 


since the photon (photino) part obviously gives no contribution to this anticommutator at the 
classical level. However, the gauge held contribution emerges at the one-loop level due to the 
quantum anomaly. This is most readily seen by using the Pauli-Villars regularization (Fig.l). 
If we add the Pauli-Villars helds to render the theory ultraviolet- hnite, all commutators 
do reduce to the naive tree-level expressions, with the regulator term explicitly added. In 
particular, in SQED 

{Qa Qp} = I d'xV + Mo{R:r;)] , (29) 


where is the lowest component of the regulator chiral superheld. It enters in the La- 
grangian in the same way as S, T, with the substitution m —>■ Mq (Mq —>■ cxd) and the 
opposite metric, i.e. the sign of the regulator loop is opposite to the normal one. 


X 


R 


* 


R 


* 


—X—^ 

\|/>^ \\f^ 

’s \ 


X 


FIG. 1. The anomalous part in the central charge. 


Integrating out the regulator held in the gluino background (this amounts to calculating 
the diagram depicted on Fig. |l]) we obtain: 


{Qa Qib} = (-4i)(a)„/3 








AA 


(30) 
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The right-hand side can be written as the lowest component of the superfield in the operator 
equation, 


{Qa Q 13 } = (-4i)(a)„/3 




mST + 




WW 


Jd=0 


(31) 


where Wa is the chiral superfield generalizing the photon held strength tensor. The expres¬ 
sion in the square brackets is nothing else but the Konishi anomaly [0 (see appendix A and 
Eq ■ (^T|) ) ■ As was expected, the central extension is determined by the quantum anomaly of 
the theory. Since it is a full superderivative, 

(Se^S + Te-^T) , 

8 '' 


the central charge vanishes in SQED, and the derivation of Eq. (^Tj) is meaningful only to the 
extent that it graphically illustrates how the anomaly (the term oc W^) can appear in Tap- 
SQED is clearly uninteresting in this context and is studied only as a technical exercise. 
Although the central charge is identically zero in this case, we still hnd it instructive to 
examine the formal operator appearing on the right-hand side of Eq. m under the given 
dehnition of the supercharges, see Eq. (P7|) . 

Let us proceed now to much more interesting case of the non-Abelian theories. In the non- 
Abelian gauge theories, the direct calculation of the central extension is more complicated, 
for it is necessary to include the contribution of the gauge helds and, as always, one has 
to hnd a way to regularize both, the ultraviolet and infrared divergences in a gauge and 
superinvariant way. It is easier to express the anomaly in {QaQg} in an indirect way by 
exploiting the anomaly relations already established in the theory. 

For dehniteness let us consider SU{Nc) SQCD with Nf havors. Each havor consists of 
two subhavors, Qa and Qa where a is a SU(N) index. The mass term is assumed to be 
diagonal. 




(32) 


where the sum runs over all havors. Such a diagonalization is always possible. In the general 
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case more generic terms in the superpotential are allowed. Our final result will be valid for 
a generic superpotential. 

The set of “geometric” anomalies of SQCD includes that in the Rq current, supercur¬ 
rent and the trace of the energy-momentum tensor. In the superfield notation this set is 
represented by the anomaly in the supermultiplet of currents |§,§], 


n" 7 ■ = -D 

^ ^ OlOL 


1 

3' 




dW 

dQi 


^T{G)-Y.iT{R, 

Ibvr^ 




This is a general expression, valid for any gauge theory, with the arbitrary gauge group and 
superpotential. In SQCD we must use superpotential (|3^ , while T{G) = and T{Ri) = 1 
for each flavor. Equation (|33|) is valid even in the presence of the Yukawa (trilinear) couplings 
in the superpotential. 

The lowest component of the current supermultiplet Jaa is the Rq current, while the 9 
and 9 components of the supermultiplet are related to the supercurrent Jpaa, see Eq. (i- 
Now, we can obtain the central extension by investigating the commutator of the supercharge 
Qa with the supercurrent. 




(34) 


With the aid of Eq.(^, the supercurrent can be expressed in terms of the 9 component of 
the current supermultiplet. 





{'A'i'} 


0/3 


(35) 


Then we have 




(36) 


The 9 component of the anticommutators on the right-hand side can be rewritten as the 
lowest component of the super derivative of the current. 


i{QaJ'yy}0l3 SSjaD 


"77 


0=0 


(37) 


(33) 
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The last term being antisymmetric with respect to the indices a (3 does not contribute to 
the anticommutator of the supercharges and we drop it, while for the first term we obtain, 




0=0 


(38) 


Assembling everything together we get 


[QaQi,} = - 2 (< t “)2 / 


(39) 


Substituting the anomaly equation (^) into the superderivative of the current we get 


—4i 




3W-^g. 


dW 


dQi 




0=0 

0=0 


(40) 


This is our master formula. In deriving it we took into account the fact that the term 
with time derivative is proportional to and cancels out after the symmetrization over 
the indices a, jS, while those with spatial derivatives are proportional to the matrix (a) 





Equation (^0[) is a general, multiloop expression for the central extension of the = 1 
superalgebra in the gauge theories with a superpotential. If W 7 ^ 0, as was mentioned 
above, the superalgebra acquires the central extension already at the tree level. Note that 
the classical part (one containing the superpotential and its derivatives) in Eq. (|40D is in 
agreement with that obtained for the Wess-Zumino model (Eq. ([T^) ). As was expected, the 
result is the integral over a total derivative. It can only be supported by extended field 
configurations. Moreover, Eq.(^) contains total super derivatives that do not contribute to 
the central charge (i.e. to the matrix elements in the supersymmetric vacua). One such total 
derivative is explicitly present in Eq. (^OD , while the other can be seen by invoking again the 
Konishi anomaly 0, 


i674+e''Q.) 


1 dW 


J2iT{Ri) 2 

167r2 


(41) 
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Subtracting the Konishi anomaly from the right-hand side of the Eq.(^0D we obtain 


{QaQp} = (-40(<?)a/3 / 


16ir2 


-|- t.s.d. 


0=0 

0=0 


(42) 


where t.s.d. denotes total super derivatives. The expectation value of the operator appearing 
on the right-hand side over any supersymmetric state is renormalization-group-invariant, 


to all orders. This was established in Ref. ||^, see Sect. 4. This fact can serve as an 
independent prove that the central extension of iV = 1 superalgebra derived above for 
supersymmetric gauge theories is exact to all orders. 

We can now compare this result with the independent calculation of the central extension 
in the Abelian case (Eq.(^lD). Indeed, for the Abelian gauge group and superpotential (|3^ 
Eq.(^^ reduces again to the Konishi anomaly, in accordance with Eq.(^). Unlike SQED, 
Eq.(0^ is non-trivial. The non-Abelian gauge theories generically do possess discrete sets 


of vacua, giving rise to the possibility of the domain wall solution [|,^]. Correspondingly, 
the central charge need not vanish. The non-vanishing central charge in SQCD will be 
transparently seen from the consideration below. 

It is quite natural, that the central extension is expressed in terms of anomalies already 
existing in the theory. If it were not the case, it would imply the existence of a new anomaly 
of a geometric nature, not known so far, that appears only in the anticommutator of the 
supercharges. No such new anomaly exists. 

To reveal in a more graphic form the relation between the “old” anomalies known in 
the gauge theories, and the anomalous term in the central charge under consideration, we 
suggest a slightly different derivation based on the so called Veneziano-Yankielowicz effective 
Lagrangian p3|,^. The advantage of this derivation is that it translates the (anomalous) 
loop effects in a straightforward tree calculation similar to that performed in the Wess- 
Zumino model. 

From dimensional counting and the analysis of the R charges, it is clear that the most 
general expression for the central extension oi N = 1 superalgebra has the form. 


{QaQ/?} = (-4*)(a)a/3 / d^xV W{Q,Q)+CTtW^^ 


0=0 


(43) 


16 











where C is a numerical constant. Dimensional arguments and R symmetry tell us that the 
anomalous part can be proportional only to the operator TtW^. There is no other gauge 
invariant operator in the theory, built from the gluon (super) helds, with the appropriate 
dimension and R charge. Our task is determining the numerical value of the constant C. 

The anomaly content of the theory is explicit in the Veneziano-Yankielowicz effective 
Lagrangians 

If we introduce the (composite) superhelds 


5 = 


327r2 


TrW^ , 




(44) 


the effective Lagrangian incorporating all standard anomalies takes the form 

1 ^ „N 1/3 1 


d^dS In 


c = -Jd^e(s sy +^J (ttM m) + 

^iVc-AT^detM I f ,2 nrr. r n rr ^ 

- + - y d^OTr {mM) + H.C. 


^Nc-N 


(45) 


up to corrections discussed in Ref. [|^. These corrections are irrelevant for our present 
purposes. The scale factor A is put to unity, so that all dimensionfull quantities are measured 
in units of A. The specihc choice of the kinetic term in the expression above is unimportant. 
What is important is the normalization of the various terms in the superpotential. We pause 
here to make a few comments on how this normalization is hxed. First, in the massless limit. 


SQCD with Nf flavors possesses an anomaly-free R symmetry |]25| , 


0 ^ Se” , W W'e” , M‘ ^ 


(46) 


and, additionally, SU{Nf)L x SU^Nf)^ global invariance. The latter requires the effective 
Lagrangian to depend on detM, while the former unambiguously determines the argument 
of the logarithm in Eq. (|45|) . The coefficient of the hrst term in the square brackets can be 
found by inspecting, for instance, the conformal anomaly. Under the transformations 


X 


x{l + 5), G“^^G“,(l-25) ... 


fll' \ 


(47) 
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the action of the fundamental theory 




X. 


(48) 


The effective action (|^ ) has the same variation. Since it is supersymmetric, all other 
“geometric” anomalies are automatically reproduced. Finally, the coefficient of the last 
term in the square brackets is adjusted in such a way that at the minimum 




The latter expression is a consequence of the Konishi anomaly 0, 
io" (Qe''Q + Qe-''Q) = mQQ + 


(49) 


(50) 


The superpotential in Eq.(^5D has minima of two types. First, the chirally symmetric min¬ 


imum at S'* = 0, observed in Ref. [|1^. (All the quantities at the minima are marked by 
the asterisk.) The corresponding value of the superpotential is >V*i = 0. Second, it has Nc 
“standard” chirally asymmetric minima at 


= detM* = 


The corresponding value of the superpotential is 


3\^/ 


mi... niNf 


(51) 


= 


2N^ 


'-S, 


(52) 


Now, the Veneziano-Yankielowicz effective Lagrangian is nothing but a generalized Wess- 
Zumino model. The central charge is immediately obtained from the tree-level calculation, 
and is given by Eq. ([T7|) . Substituting Eq. (|52D we hnd 


S = 2 [>V*2 - >V*i] = 


ANr 




(53) 


The latter result is compatible with the general expression 

Nf-N, 


C = 


Ibvr^ 


at 


(54) 
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provided one takes into account the Konishi anomaly (|^ . This value of constant C is in 
agreement with Eg. (^2D . 

Although, at hrst sight, it may seem that the derivation above relies on the presence of 
the additional, chirally symmetric vacuum at 5* = 0 (which was introduced in Ref. [|^; its 
existence is still to be conhrmed), actually the calculation with minor changes stays intact if 
we consider the domain wall, that interpolates between two “standard” chirally asymmetric 
vacua. The hnal result for the constant C is unchanged. It depends only on the general 
symmetry structure of the Veneziano-Yankielowicz effective Lagrangian. 

Let us return to SU(Yc) SQCD with Nj flavors and ask the question why Nf = Nc is 
special so that at Nf < Nc the coefficient C in front of the anomaly does not vanish, while 
at Nf = Nc it does. At Nf < Nc a superpotential for the (composite) held M is generated 


non-perturbatively even in the absence of the tree-level superpotential W [25 


Wn.p. ~ 


' Jy3Nc—Nf\ Nc-Nf 


detM 


(55) 


In our language this non-perturbatively generated superpotential is entirely due to the sec¬ 
ond, anomalous term in {QaQf^}- At Nf = Nc the conserved i?-charge argument forbids 
non-perturbative superpotential [^. The vanishing of C at Nf = Nc is in one-to-one 
correspondence with the absence of any function of detM that would have appropriate R 
charge. 

It is clear, that two derivations of the central extension presented above are not inde¬ 
pendent. They both relate the anomaly in the anticommutator of supercharges to other 
anomalies of the theory. A direct computation of the anticommutator, say, by using the 
Bjorken-Johnson-Low limit 0. is more complicated due to complications with the gauge 
invariant supersymmetric regularization. Thus, the full diagrammatic calculation of the 
central extension is still missing. We have no doubts, however, that when it is carried out, 
our result (^OD will be reproduced. 
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III. BPS-SATURATED SOLUTIONS IN THE GENERALIZED WESS-ZUMINO 


MODEL 

The generalized Wess-Zumino model is the theory of the type (Q) with an extended set 
of the interacting chiral superfields i = 1 , 2 ,..., and the superpotential W(<hi), which 
will be assumed to be analytic function of the superfields and various parameters (coupling 
constants, mass terms and so on), if not stated to the contrary, see Eq. below. Typically, 
one can think of W($i) as of a polynomial. In certain instances of physical interest, however, 
one has to work with non-polynomial W(<hj), see e.g. Ref. [^. 

In the previous section we demonstrated that in many instances the extended classical 
field configurations lead to the existence of the central extension of iV = 1 SUSY superalgebra 
(!)■ 

The central extension turns out to be very important for the analysis of the theory in 
the sector with the domain walls or other topological defects, e.g. strings. As was already 
mentioned, the domain walls ruin the invariance of the theory with respect to the Poincare 
group, and, with necessity, lead to (at least partial) supersymmetry breaking due to the 
algebra (0. The question whether or not supersymmetry is broken completely was shown 
1]^ to be related to the appearance of the central extension. If the central charge vanishes the 
domain wall necessary completely breaks N = 1 supersymmetry. With a non-zero central 
charge two out of four supersymmetry generators may remain unbroken. 

A. Preservation of 1/2 of SUSY and BPS saturation 

In Sect. II it was demonstrated that if the domain wall under consideration is BPS- 
saturated, the domain wall conhguration preserves a part of supersymmetry. In the begin¬ 
ning of this section we will demonstrate that the inverse assertion is also true: 

Theorem 2. If a part of supersymmetry is preserved for the given domain wall, it is 
BPS-saturated, and its energy density £ = |S|. 
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Let us consider the generalized Wess-Zumino model, with a set of interacting chiral 
helds. Although we use a specihc model, the results will be quite general, and, in particular, 
applicable to the gauge theories. The Lagrangian of the model has the form 

C = d^ed^x<!>A + [I J d^9d^xW{<!>^) + H.C.^ . (56) 

We assume that the superpotential W is chosen in such a way, that the scalar potential has 
more than one minimum 


dW 

d(j)i 


0 


(57) 


at 


{0i} ) {0i*}2 ) • • • 

where, as usual, the asterisk marks the values of the scalar helds 0* at the minimum, and 
the subscript 1,2, ... indicates that the set of minima includes several (isolated) points. 
To ease the notation, in what follows, we will omit the index i, referring to the scalar helds 
at the minima number 1 ,2 ,... as , 0*2 and so on. 

The boundary conditions are 

{ 0*1 at z ^ —oo 

(58) 

0*2 at z ^ +00 

The energy of the static held conhguration with the given boundary condition (the 
domain wall) is 


S 



dz 




where V is the scalar potential 


v = E 

i 



(59) 


(60) 


Since the energy is the sum of two non-negative terms, it can always be minimized. The 
held conhguration with the boundary conditions (|5BD minimizing the energy £ is the domain 
wall solution. 
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(61) 


= 


d 


d<Pt 


E 


dW{<p) 


d(t)j 


The domain wall solution, which we have chosen to be in the xy plane, always exists since 




(62) 


It is convenient to denote the derivative with respect to z as a dot over the corresponding 
letter. Moreover, we will arrange a set of 0’s as a “vector”, 


(63) 


\---J 


Then Eq.(|61|) takes the form 


= -VW 


(64) 


where 


V = 


d d 


’ 502 + ’ ■ ■ ■ 


the dot over letters denotes the derivative with respect to and 


w = -E 


dW 


dcpi 


(65) 


( 66 ) 


Equation (|6^) is nothing but a complexihed version of the Newtonian equation for the 
mechanical motion in the potential U. For purely real solutions Eq.(^) reduces just to the 
regular Newton equation. The prohle of U can be viewed as a mountain ridge with the 
summits at 0*i , 0*2 etc. The height of all summits is the same. The analog motion starts 
at one summit and ends at another. 

So far the only new element introduced in the standard theory of the domain walls by 
supersymmetry is complexihcation. A more novel feature to which we proceed now is the 
existence of the special solutions satisfying the hrst order differential equations, rather then 
the second order equation (|^ . 
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In the spinorial notation the supersymmetry transformation has the form 


= V2^F+ - y2i(a^0+)(a^)““e„ . (67) 

Assume that (j)'^ depends only on z (i.e. 0+ is independent oft,x,y, the (x|/)-plane oriented 
domain wall), and 

Tie = ie*“e (68) 


where e and e are two-component columns 


e = 




V" y 




e = 


^ 2 * 


(69) 


V' / 


and the factor +i is singled out for convenience. (The matrix ri appears in Eq.(^) since 
both e’s in Eq. (|68D have upper indices; therefore, in Eq. (|67|) we must substitute 

e/3 = —'*(T2)/3ae“ , 


see Appendix A, Eqs. (|Allj) , (|A12|) ). The solution of Eq. (| 68 |) is 


eo = 


A i 




e*2. 


(70) 


Then = 0, provided that 


= -F+e-*“ = 


dW 
d(j) ^ 


(71) 


If equation (^) is satished, two out of four supercharges will annihilate the domain wall, the 
corresponding part of supersymmetry will be preserved More exactly, Eq. (|7TD must 

hold for all helds 0 / from the set, with one and the same value of the phase a. 

Under the standard choice of the superpotential 


>V*2 = >V((/.*2) = >v [C(^ = +oo)] = K ( 72 ) 

where k is real and positive (obviously 2k = |S|). This choice is convenient since under this 
choice we can £x the phase a in Eq.(^) at zero. 
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In the spirit of our mechanical analogy we can rewrite Eq. (O in the following way 


(f=VW. 


(73) 


The “gradient” is dehned in Eq.(|65D. This is nothing but a (complexihed) equation of 
motion of the high-viscosity fluid whose inertia can be neglected. We will refer to Eqs.([73D 
as streamline or cree/sequations. The creek flow starts at one saddle point of VV (i.e. >V(0*i)) 
and ends at another saddle point (VV( 0 * 2 ))- 

It is quite clear that if Eq. (|73D is satished the second-order relation ( 0 ) is satished too. 
Q Indeed, start from Eq.(|73D and differentiate both the left- and right-hand sides. We then 
get 


V _ ^ d^W dw 
dz ^ 92; \ d(j)f J j 90+90+ ^ ^ 90+90^ 90j 

which identically coincides with Eq. (|^ . The opposite is obviously not always true: not 
every solution of Eq.(l^) is simultaneously the solution of the creek equation ([7l|) . 

Now we are hnally able to prove Theorem 2 formulated above. Assume that we deal with 
a solution of the creek equations which then automatically preserves one half of SUSY. The 
energy of the held conhguration is (see Eq. (|5^) ) 


/ +00 

dz 

-OO 


E' 


9W 


2A (W*2 - W*i) = A |S|. 


(75) 


The fact that 8/A = |S| means, by dehnition, that the wall is BPS-saturated, q.e.d. 

For completeness we give here an alternative proof (Refs, 
plementary to that given in Sect. II.B. Consider 


16| ) of Theorem 1 , com- 



(76) 


The integral o is obviously positive-dehnite for all 0*. At the same time, it can be rewritten 
as 


* Analogous equations in non-super symmetric models were considered in Ref. |]T] 


24 















■ + 00 


(77) 



dw a>v \ 

d(pt dcpi) 



> 0 . 


Minimizing the first term by varying (f)i we get 


£ > |S 


and the equality demands the domain wall function to satisfy 



Thus, this is necessary and sufficient condition. 


B. Properties of the BPS-saturated walls 


Within our convention regarding the superpotential we can prove a stronger statement. 
Theorem 3. Let the BPS-saturated domain wall exist. Then W [0^(2^)] is a real positive 
number for all provided our standard convention on the superpotential is adopted. 

This relation is a direct consequence of Eq. (O)- Indeed, let us multiply both sides by 0 
and integrate over z in the interval from —oo to the current value of 2 ;. Then we get 


{>v[0rw]->v.i}= r d^i0p>o. 

J —00 


(78) 


This proves the theorem, since according to our standard choice of the superpotential >V*i = 
0. A particular case of this result was used previously p2[. 


The superheld describing the domain wall so far has only bosonic components, 
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+ 0 = d)^ - 9 


dcp 


(79) 


If we add the fermion zero modes, generated by those supertransformations that act on the 
wall non-trivially, we will get the wall superheld, which has all three components, the lowest, 
the middle, and the F component. In the case of the BPS saturation the wall superheld 
satishes the full supersymmetric equation of motion. 
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a>v($) 

d^i 


( 80 ) 


for further details see Sect. IV.A. 

It is obvious that if the central charge E = 0 the corresponding domain wall breaks 
supersymmetry completely |]^. From Theorem 1 it is clear that the BPS-saturated solutions 
preserving a part of SUSY can exist only provided E 7 ^ 0. Unfortunately, the non- vanishing 
central charge is not sufficient to ensure the existence of the BPS-saturated domain walls. 

Let us examine the example suggested in Ref. [|l|: the generalized Wess-Zumino model 
with two chiral superhelds, X and <F, and the superpotential 


W ^x(i? - 


(81) 


(cf. Eq.(^). This superpotential has two minima: 


= -^,X = 0 ^ and {(() = 


/i 


= 0 


(82) 


^/A’" j r %/A’ 

where y is the lowest component of the superheld X. The corresponding values of the 
superpotential at the minima vanish. Hence, the central charge E = 0 in the model (pTD, 
and the creek equations have no solution. At the same time, the classical equations of motion 
( 1 ^) for the held 0 identically coincide with that of the model (^ provided we put X = 0 . 
The solution is of course the same, as in the model (0), with the energy density 

8 / 1 ^ 


£ = 


3v/A 


(83) 


(See Sect. 4 for further details; it is assumed for the time being that fi and A are real positive 
parameters.) As we know, the solution of Eq.(| 6 T|) always exists. The fact that £ > 0 while 
E = 0 shows that the domain wall is not BPS-saturated, and, correspondingly, SUSY is 
completely broken. 

Now, let us deform the model in such a way as to generate a non-vanishing central charge. 


>V(a) = X - A4>2) + . 


(84) 
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The values of the scalar fields at the minima become 


and 


At the minima 



a /i I 

^71J 


a I 


^ -f 7 a • ^ ^iTa 

The central charge 

E(a) = 2AW(a). = # 0, (86) 

(the parameter a is assumed to be real and positive ). The domain wall solution of the model 
(0) will also be deformed. It is clear, however, that as long as a/A -C 1 the deformation is 
small, and the energy density e{a) is close to Eq.(|^). Thus, for sufficiently small a/X 


e{a) 3> |S(a)| , 


(87) 


and the domain wall can not be BPS-saturated. 

As a matter of fact, it is not difficult to prove a stronger statement. Assume that at 
some positive value of a the domain wall in the model at hand is BPS-saturated. Then it 
satisfies the creek equations 


X+ = /i^ - , 0 + = - 2Xx4> ■ 


Multiply the first equation by <p and integrate over z from —cx) to -|-cx). 


h = 


r+oo 


dz(l)x^ = 


r-\-OQ 


dz 


7. 

dz 


/iV- 


A 


3^ 


The integral 


h 


/oo {i - x) (i* - X*) = 


( 88 ) 


(89) 


( 90 ) 


27 








is obviously non-negative. The first integral on the right-hand side 


-l-CXD 


+ XX'^) dz 

If the domain wall is PBS-saturated then 


X r + OO 
2 7-00 


+ XX + 


v) dz 


e 

2 ‘ 


(91) 


£ = |S|. 


(92) 


Hence in this case we must have 


or 


|E| > 2/i 


2a ^8 /i^ 


(93) 


(94) 


3 SyX' 

We conclude that the domain wall at hand can not be BPS- saturated at a/A < 4. If there 
is any chance to make it BPS-saturated a/\ must be > 4. 

Actually, there is no chance. The theory under consideration has two fermion helds, V'# 
and 'ipx- Consider the fermion mass matrix /d^id^j in either of two vacua (pS]). After 
diagonalization of the mass matrix we hnd two independent diagonal linear combinations of 
"0$ and 'ipx (they are the same in both vacua). The sign of the mass term changes in passing 


from one vacuum to another. According to the index theorem the Dirac equation for 
each of the diagonal combination will have a zero mode. At the same time only one linear 
combination of the bosonic helds has a zero mode. Hence, SUSY is broken completely on 
the domain wall for any a. The absence of the BPS-saturated solutions for any a for real 
values of the helds is readily seen from the consideration of the prohle of W. 

Conjecture. If in the generalized Wess-Zumino model with the polynomial potential 
the domain wall is BPS-saturated at some values of parameters in the superpotential, it 
continues to be BPS-saturated for all non-singular deformations of parameters as long as 
the central charge remains hnite (neither inhnity, nor zero), and 

1 


det 
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does not cross zero. If at some values of the parameters the domain wall is not BPS- 
saturated, it can not become BPS-saturated under arbitrary non-singular deformations of 
the parameters. (Adding to the superpotential terms of higher powers in the helds which 
were originally absent is a singular deformation). The latter statement is valid provided 
there is no level crossing, i.e. the energy densities of two domain walls, £i(A) and £ 2 (A), 
do not collide at some A (here A is a generic notation for the set of parameters in the 
superpotential). 

The simplest example of the BPS-saturated domain wall is provided by the minimal 
Wess-Zumino model, see |l|] and Sect. 4. 

Unfortunately, we were unable to hnd a sufficient condition for the existence of the 


BPS-saturated domain walls in the general case. In Ref. the BPS-saturated walls were 
found in the Veneziano-Yankielowicz Lagrangian which gives an effective description of the 
supersymmetric Yang- Mills theories (with or without matter). 

If the central charge is non-vanishing, the parameters of the superpotential are purely 
real, and one decides to limit one’s searches for the BPS-saturated domain walls by purely 
real solutions, one can take advantage of the fact that Eq.(|^) is analogous to the one 
describing the motion of the high-viscosity liquid on the multidimensional profile given 
by the superpotential W. This mechanical analogy and our intuition in solving mechanical 
problems, helps in establishing the condition for the existence of the BPS-saturated solutions 
in this case. Making use of the mechanical analogy, we can infer that if the superpotential 
has a set of extrema, and at least one of them is either maximum or minimum, then there 
always exists a BPS-saturated domain wall solution, and 1/2 of SUSY is preserved. (If all 
extremum points are saddle points, the BPS-saturated solution may not exist.) 

It is instructive to present here the corresponding example - a model in which the 
BPS-saturated wall obviously exists for real values of the parameters in the superpotential, 
and it remains to be BPS-saturated when one analytically continues the parameters of the 
superpotential in the complex plane provided that the central charge stays finite (does not 
go to zero or infinity). The model we keep in mind has the following superpotential 
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(95) 


W(*, X) = -oA' + ^4'^ + |(A - 4)“, 

Zi o 

where a, (3 and 7 are parameters. To begin with, we choose all three parameters real and 
positive. Then the solution of the creek equations can be sought for in the class of the real 
functions. By inspecting the prohle of W it immediately becomes clear, that this prohle 
has one maximum and one saddle point, and the creek equations do have a solution. The 
analytic form of the solution is unknown, but it exists, for sure. Moreover, the extrema of 
W are 

a a la 

<P*i = ^ y X*i = ^ 


P 


P 


7 


(maximum) and 


a a a 

0*2 "7 ) X *2 "7 \ 

p p \ 1 


(the saddle point). The central charge 


a 


E = 2AW, = -aj- 
3 V 7 


(96) 


is non-singular and non-zero provided a and 7 7 ^ 0, cx). The determinant of the matrix of 
the second derivatives of the superpotential at the extrema is 

207^^, and - 207 ^^, 

respectively. We can now move the parameters into the complex plane keeping 0 7 ^ 0; the 
BPS-saturated wall will persist. The wall energy density will be determined by the absolute 
value of ( 1 ^. 

In summary, the existence of the central extension with a non-vanishing central charge 
is a necessary condition for domain wall to be BPS-saturated |l|. The absence of the 
central extension (i.e. vanishing central charge) necessarily leads to the complete breaking 
of supersymmetry, while in the presence of central extension supersymmetry may or may 
not be broken completely. 
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C. Three or more vacua 


So far we essentially disregarded the possibility that one and the same model can have 
more than two vacuum states, and, correspondingly, a variety of distinct domain walls. 

Let us start, however, with a remark referring to models with one pair of supersymmetric 
vacua. One and the same model can have simultaneously two types of the domain walls 
interpolating between the same vacua. The picture becomes transparent if one limits oneself 
to superpotentials with real parameters and real solutions of the creek equations. 

The hrst type of domain walls appears as the solution of the classical (second-order) 
equations of motion (^) in the bosonic sector. These solutions do not satisfy the creek 
equations. It is clear, that supersymmetry is not preserved on such walls. Solutions of the 


second type, are those of the hrst order differential equations (|80D, and they are partially 
supersymmetric by construction. Classically, both types of solutions may coexist peacefully, 
and be stable. The energy density in the hrst class is higher than in the second, and 
they become unstable quantum-mechanically. The “false” domain wall will transform into 
the BPS-saturated through the formation of “bubbles” on the wall which will grow in an 
explosive manner. 

If there are more than two vacua, the picture may become even more versatile. Some of 
the vacua may be connected by BPS-saturated walls, others by non-saturated walls. Let us 
assume for simplicity that the model at hand has three vacuum states, and all are connected 
by the BPS-saturated walls. We would like to consider a superposition of two walls - a 
held conhguration that hrst interpolates between the hrst and the second vacuum, and then 
interpolates between the second and the third. If these two walls are separated by a hnite 
distance in the 2 ; direction, the corresponding held conhguration is not the exact solution (it 
only approaches it when the separation tends to inhnity). The behavior of the superimposed 
walls crucially depends on the corresponding central charges (which are generically complex). 
Let us denote by Sjj the central charge for the wall conhguration interpolating between the 
Tth and j-th vacua. If all Sjj lie on one line in the complex plane and are ordered along 
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this line then, obviously, the superposition of walls 12 and 23 has higher energy than the 
BPS-saturated wall 13 would have. However, the infinitely distant walls 12 and 23 will have 
exactly the same energy as the BPS-saturated wall 13. This means, that there is a repulsion 
between walls 12 and 23 which will eventually push them out at infinite separation. If 
there is an independent BPS-saturated interpolation between 1 and 3, these two types of 
conhgurations will coexist. 

The above case is clearly degenerate. If the central charges do not lie on one line and 
form a triangular, with the corresponding inequality, 

IH12I -l- IH23I > iHial 


the infinitely distant walls 12 and 23 will have larger energy as the BPS-saturated wall 13. 
This will lead to the wall attraction at large distances. The walls 12 and 23 will collapse, 
forming a lower e state, the wall 13. 

It is interesting that many elements of this rich picture are actually realized in supersym¬ 
metric SU{N) Yang-Mills theory |^. In this model there exist N domain walls connecting 
N chirally asymmetric vacua with the chirally symmetric one |]^ at (AA) = 0. Let us call 
them primary. These walls are trivially BPS-saturated. In the SU{2) model there are three 
vacua, and all central charges are real. Hence, if we superimpose two primary walls, we 
will get a repulsion at large distances. For SU{3) and higher groups SU{N) there are also 
N{N — l)/2 walls connecting directly various chirally asymmetric vacua, which also seem to 
be BPS-saturated The primary walls experience strong attraction at large distances; 
the walls connecting the chirally asymmetric vacua can be viewed as bound states of the 
primary walls. 


D. BPS-saturated strings 

The very existence of the walls is due to the fact that our space is not compact - the 
dynamical fields at all infinitely distant points need not have one and the same value. We 
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can use the very same idea to consider other topological defects, in the spirit of what was 
done in Ref. []^. Our consideration here is conceptually close to that given in Sects. 3 and 
5 of Ref. 1^, but we will focus on those defects which are BPS-saturated, i.e. 1/2 of the 
original SUSY is preserved in the corresponding background held. The conhgurations to be 
discussed are string-like solutions, static and z independent. They spontaneously break two 
out of four Lorentz translations, so that the resulting dynamical system has H = 1 -|- 1. 

Let us orient the string along the z direction, so that all bosonic helds in the solution we 
are going to build will depend only on x and y. Introduce 

(^ = X + iy and ( = x — iy . 


The central extension of the superalgebra, see Eq. (0, now takes the form 

1 


{Q^ Qf] = 4 L y dC 5(1 - 


w|4t(C.C)] 


1^ , aw|$t(c,OI 

3 T ‘ 94. 


'K 2(1 + 


w|4t(C.O] 


i^^- 9W|4t(C,C)| 
3 T ‘ 94. 


(97) 


e=o 


where L is the length of the string, and the integrals run over the large circle in the ( plane. 
The expression in the braces is the central charge. This expression assumes that the helds 
<hj are non-singular at hnite distances, and the superpotential W is non-singular at hnite 
values of the helds. Introduction of the gauge helds results in obvious modihcations, with 
the anomaly term appearing on the right-hand side. 

What supercharges are conserved on the BPS-saturated strings? To answer the question 
we must go back to the general supertransformation law (|67|). If the parameter of the 
supertransformation is chosen as follows 


e = 




V' / 


/ 


0 




exp(z7r/4) 


(98) 


then S'lp = 0 provided 


1 dW 

2 debt 


(99) 
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Here 



and 



We could have introduced an arbitrary phase in Eq. much in the same way as in 


the discussion following Eq. (^). This phase can always be “unwind” by an appropriate 


phase rotation of the superpotentiah Needless to say that interchanging the lower and upper 
components of e we would arrive at the first-order equations defining the “anti-string”. 

We will continue to call (p^) the creek equations, although the “time” variable is now 
complex. 

It is seen that the creek equations (p9|) corresponding to BPS saturation in the string 
case, are modified in a minimal way compared to the domain wall problem. As in the latter, 
there are several alternative lines of reasoning allowing one to obtain the creek equations. 
Say, we could have started from the analog of Eq. (|76D and consider 



( 100 ) 


This integral is obviously positive-definite; it becomes zero only if each of the brackets 
vanishes. This is the condition of the BPS saturation. If Eq. (p9| ) holds, then the standard 
equations of motion obviously hold too. Indeed, 



( 101 ) 


Up to a total derivative / reduces to the string tension a minus the central charge. If 
the latter is non-zero, the hrst-order equations may have a solution. Then the string 
tension will be exactly equal to the central charge. 

The statement above requires an immediate reservation: I reduces to the string tension 
minus the central charge provided we discard total derivatives. Presumably, this can be 
achieved by adding certain total derivatives directly in the Lagrangian. The issue is not 
completely clear at the moment, further investigation is needed. 
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Let us consider a particular example, the model with the superpotential 


W = X - 

72 

For simplicity the mass parameter in front of X is put to unity. In the conventional sector, 
without topological defects, this model has a runaway vacuum, 

v = i|xY + |i-72^<h|7 



corresponding to df ^ 0, <F —>■ oo, much in the same way as in SU{2) SQCD with one flavor 


and the vanishing matter mass term. It is easy to see that the creek equations (|^) have 
the solution 


2^ = C(i + C0-7 


1 

72 


C 1 + (1 + CO 



(103) 


This solution has a constant volume energy density due to the kinetic term of the $ held, 
as in Sect. 5 of Ref. |]l|. This depends on the dehnition (i.e. on whether or not we add total 
derivatives). Discarding this constant volume energy density we observe a typical string-like 
distribution of energy, with a hnite string tension. The central charge is non-vanishing since 
at large |C| the linear part of the superpotential W behaves as 1/7 (The cubic part can be 
ignored as is evident from Eq. (^)). 


IV. MINIMAL WESS-ZUMINO MODEL 

In this section we carry out the construction of supersymmetric theory in (2 -|- 1) di¬ 
mensions, obtained by dimensional reduction of the minimal Wess-Zumino model in (3-1-1) 
dimensions on the domain wall solution. We construct an explicitly supersymmetric La- 
grangian with the inhnite set of interacting superhelds. Each superheld realizes an irre¬ 
ducible representation of iV = 1 supersymmetry in (2 -|- 1) dimensions. We, then, compute 
the quantum corrections to the shape of the domain wall. In Sect. 5 we will use this informa¬ 
tion to obtain the hrst subleading correction in the asymptotic behavior of the multiparticle 
production amplitudes at threshold when the number of particles tends to inhnity. 
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A. Classical solution 


The Lagrangian for the minimal Wess-Zumino model was written in Eq. (^. For our 
present purposes it is more convenient to define the parameters in the superpotential as 
follows: 


W(4) = ^ 4 - jit 


(104) 


where 


* = (> + + eoF 


(106) 


is the only chiral superfield of the minimal model. It consists of one complex scalar field 0, 


one chiral fermion -0, and an auxiliary field F. 

Although at first sight it might seem that the choice of the superpotential is not general, 
actually it is not difficult to see that any renormalizable model with one chiral superfield 
can be reduced to Eq. ( |104|) by exploiting the freedom we have in defining the superfield. 
After $ is appropriately shifted and rotated we arrive to Eq. (|104j) , with the real values of 
parameters, no matter which cubic superpotential we started from. 

Lagrangian (H) has a discrete Z 2 symmetry group (cf. Eq.(||)), which is spontaneously 
broken, so that there exists a classical solution of equations of motion: 



4 


W 


a>V(<f>) 


(106) 


where 


$. = 0 . + eoF^ = cj)^- eeW{<pt,), 


(107) 


and F = -W = -dW/d(t)+. 

Here 0^, is the real field 0^ = 0+ 



( 108 ) 


satisfying 
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= W'(M . (109) 

OZ 

The domain wall was chosen to lie in the xy plane. 

Note that in the minimal Wess-Zumino model (when the solution is sought for in the class 
of real functions 0 ( 2 :)) it is always possible to reduce the second order differential equations 
of motion (^) to the creek equations, by exploiting the integral of motion (conservation of 
energy). Thus, the fact that the domain wall in the minimal Wess-Zumino model preserves 
1/2 of SUSY is rather trivial. The domain wall is, then, BPS-saturated and its energy 
density coincides with the topological charge 


£ = IS 


8m^ 


( 110 ) 


Let us discuss now the “other” half of supersymmetry, namely those supergenerators 
that are explicitly broken by the wall. The domain wall at hand satisfies the creek equation 
( p.09|) (the wall profile is purely real). The generic supersymmetry transformation (^) in 
this case can be identically rewritten as 


50 = -V 2 


dz 


(e + ine ), 


( 111 ) 


where in the remainder of this subsection all spinors are assumed to be taken with the upper 
indices (which will not be shown explicitly). 

Applying the broken symmetry generators to the domain wall superfield 0^ produces 
zero modes. For instance, action of the broken generator of translation along the z axis 
gives the bosonic zero mode. We are interested now in the fermion zero mode. The super¬ 
transformation generating the fermion zero mode corresponds to the parameter 


e = 


U 1 


V“V 


( 112 ) 


so that Tie = —ze, and 





(113) 
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To describe simultaneously all zero modes, fermionic and bosonic, it is convenient to 
introduce the wall superfield much in the same way as it was done for the instanton 
superheld in Ref. (see especially Sect. 3 in this paper). To this end we hrst need to 
introduce the collective coordinates. The bosonic collective coordinate is zq, the position of 
the center of the wall. In the consideration above it was assumed that the wall center is at 
the origin. Clearly, this is not the generic case; the wall center can be at any point on the z 
axis. To get the wall superheld it is necessary to restore zq- Now, we additionally introduce 
the fermion collective coordinate with the transformation property 

600 = e + *Tie . (114) 


Although this collective coordinate is formally written as a two-component spinor, actually 
it represents one complex Grassmann number. 


00 = 


( \ 

V 

\'V 


(115) 


and St] = + ie^. Needless to say that there exists also a conjugated collective coordinate. 


00 . 

The transformation properties of 0o are rather peculiar. Actually, for building the wall 
superheld we will only need to consider the diherences zr — Zq and 0 — 00- (Naturally, 
in the conjugated superheld one deals with zl — zo and 0 — 0o.) 

Then, 


6(0- 6>o) = -ine, 


(116) 


and 


S{zr - Zo) = 2i[eT3{0 - 6 * 0 )]. 


It is easy to check that the invariant interval has the form 


(A^)in^ = (zr - Zo) - i{0 - 0 o)t 2{0 - 0o) ■ 


(117) 


(118) 
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Now, if we take the domain wall solution 0+ of the pre-supersymmetry era, given in Eq. 


(p.08|), and substitute there ^ by (Az)inv, we will get the full wall superheld 


m 


^ w { zr ,9) = ytanh[m(A^)i 


(119) 


which contains the classical solution itself, and the zero modes. The fermion zero mode 
emerges when one expands (Az)inv appearing in tanh above, keeping the term linear in 9. 
The term quadratic in 9 will automatically produce the component. The 9q part of 
(Az)inv corresponds to shifting the domain wall center. The expression for is similar. 
The meaning of the parameters 9q,9q (fermionic collective coordinates) is now clear. 


B. Quantum fluctuations 


To obtain equations that govern the quantum fluctuations around the classical solution 
we rewrite Lagrangian (j^) in components and expand it in “quantum” parts of the helds 
around the classical conhguration, keeping only the quadratic in the “quantum” helds terms, 

C = J + {d^bf - {a? + 6^) [W"(0^)]^ - (a^ - 6^)>V'(0^)>V'"((/)^) 

- ^>V"(0^) ('0-0 + 00) |, (120) 

here 

a = Re(0 — 0^); b = lm(0 — 0^„). 

Potential terms are different for real and imaginary parts of the scalar held, for the 
classical background is real. This gives the following set of equations of motion: 


{d'^ + O^O\z)}a = 0, 

{d'^ + 0^0\z)}b = 0, 

0‘^{z)'^~ = 0, 

id^-+ 0\z)'^+= 0 (121) 


where we introduced a short-hand notation 
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0\z) = {d, + W'{(t)^)}, 

o\z) = {-d, + w"{K)}. 


tl)a + ^(o-^)aa'0“ 

vl>~ 

i’a - ^(cr^)a«'0“ 

L J 

5 

iV2 


g0 = (^0^3 _ ■q1^2 ^ 

& = {^o,^^,^2} (122) 

where (r*)“^ are the Pauli matrices. 

Since the classical solution 0^ depends only upon coordinate z, the 2 ; dependence of the 
helds can be factorized 


n 

=Y.^rian{z) d'+(x\ 0:°), 

n 

K^) = H a;°), 

n 

'^nix) =J2V^^n{z) (123) 

n 

where the factors a/^ or are inserted for proper normalization ( A„ and Kn are the 
eigenvalues, see Eqs. (|124|) , (|127|) below). 

Functions an{z), bn{z) satisfy the following equations: 

(D (^z'jOiYii^z'j Kj^bjii^z^^ 

0‘^{z)hn{z) = \nan{,z). (124) 

They are eigenfunctions of two different operators 

0‘^0^{z)an{z) = mlan{z), 

0^0^{z)bn{z) = mlbn{z), (125) 


where m 
Eq.dH) 


^ = XniZn, and, thus, can not be normalized to unity simultaneously. Making use of 
we easily obtain (brackets stay for the integration over z) 


^n{ bn bn ) A^( (Xn (Xn )) 


( 126 ) 
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and the normalization condition is 


1 1 
{^n^n) 5 (®n®n) T • 

'^n 


(127) 


We snbstitnte the expansion (|123|) into the Lagrangian 


, integrate over 2 ; to obtain 


three-dimensional Lagrangian for the flnctnations on the snrface of the domain wall. Each 
given mode of Eg . (11241) will correspond to a particle with the mass given by the corresponding 
eigenvalne. Since we carried ont the expansion aronnd the classical solntion only np to the 
second order in helds we obtain a free theory of two bosonic and fonr fermionic components, 


where 


C = fd^xJ2\{d^anf + {d^bnY - ml{al + bl) 

+ ]^an{id - mn)an + + rnn)l3r. 


= 


( 128 ) 


4/+ -I- *4/“ 

^ n ' ^ ^ n n ^r) v ^ ^ 

') Pn 


and d is given in Eq. (|122|) . The snmmation convention nsed thronghont the section is 


(129) 


aa = a^afs] ada = a^dja.y. 

The Lagrangian (|128| ) is obvionsly snpersymmetric. The pairs {d,a) and (fe,/9), where the 
tilded helds are some orthogonal linear combinations of the helds a and b, form two irre- 
dncible representations of iV = 1 SUSY in d = 2 -|- 1 dimensions, that consist of one real 
boson held and one real two-component spinor. 

Since SUSY is not broken completely we can try to constrnct the Lagrangian, where the 
remnant snpersymmetry is manifest not only for qnadratic hnctnations. We want to take 
into acconnt interaction between diherent modes on the wall, and make snre that it is indeed 
snpersymmetric. To this end we introdnce three-dimensional snperspace, which consists of 
two real Grassmann coordinates and nsnal (2-1-1) Minkowski space. The snperheld is 
constrncted from one real boson and two component real Majorana spinor as follows. 


$ = 0 + + eef 


(130) 
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where / is an auxiliary field. Covariant derivative is given by 



(131) 


where 7 ^ are two-by-two real Majorana matrices, 


7o = ^t- 2 , 7l = '^3 , 72 = , 


(132) 


and T are the Pauli matrices. The most general renormalizable Lagrangian that describes 
the interaction of two sets of superfields is 



where 


Kf = ii + + 97®. 


(134) 


Mass terms for fields and T*® have a different sign in accord with Eq.(|128|). In (2 -|- 1) 


dimensions the chiral transformations are absent and the sign of the fermion mass can 
not be made positive. With the aid of formulae given in Appendix B, we can cast the 
Lagrangian after integration over z, in the form (|133|) . This explicitly demonstrates, 
that the dynamics of the modes on the domain wall is supersymmetric. The absence of 
the corrections to the vacuum energy, inherent to the supersymmetric theories, ensures that 
the domain wall energy density is not renormalized, implying the domain wall to be the 
BPS-saturated state. 


C. Correction to the domain wall shape 


While the domain wall energy density is fixed by its classical value, the shape of the 
solution is subject to corrections, regardless of the presence of SUSY in Lagrangian (|133D . 
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Indeed, the “tadpole” diagram, that arises from the cubic terms of Lagrangian (|133|) , gives 


rise to the linear term in the effective Lagrangian. It shifts the equilibrium position of the 
helds and each mode thus acquires a vacuum expectation value, specihc for every given 
mode. This leads to a change in the soliton shape, 

+ (135) 

where a* ( 6 *) is an eigenfunction of Eq.( p.24|) and <pf is the vacuum expectation of 
the lowest component of the superheld Once we established the existence of 

the “tadpole” graph, it is easier to compute it directly in four dimensions. The part of 
Lagrangian with the cubic coupling has the form 


C 

1 

2 


= - y {W"'(0^)W"(0^)a (a" + b^) + 


(136) 


To compute the tadpoles we need to know the propagators of the helds entering the La¬ 
grangian ( p. 20 |) , (|136|) in the domain wall background. At this point we may notice, that 
there is a partial cancelation between the boson and fermion loops. Indeed, the following 
relation holds between the traces of the bosonic and fermionic Green’s functions, 

2W" 

Tr G^,+ii,-[x,x] = ^ ^201 = 4 >V'(0^) Ga[x,x], 

2W" 

Tr G^-^+[x,x] = ^ ^ Gb[x,x], 

Tr G^+,i,+ [x, x] = Tr [x, x] = 0, (137) 


where trace is taken over the Lorentz indices. This relation makes the cancelation trans¬ 
parent. Taking into account combinatorial factor 3 in the vertex, the total result is the 
difference between the contribution of loops with held a and b inside. 


The corresponding Green’s functions are easily found and we end up with the fol¬ 
lowing term in the ehective Lagrangian, 


^£e// ^ _ J ^4:X‘^(j)y^{z){G'^[x,x] — G^[x, x])| a(x). 


(138) 
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where is the propagator of the held a {b) in the background of the domain wall. 

The presence of the linear term in the Lagrangian shifts the equilibrium positions of the 
helds away from zero. To hnd a new vacuum conhguration we must solve the equation of 
motion for the boson held a, and, thus, obtain the corrected shape of the domain wall, 

f 1 A 

= (t)w{\ rh) + I I 0^(1 - [—(139) 

where 


A = A 1 - A^ 


Stt^ 


In(Mo) 


87rv/3 


m = m \ 1 — 




ln(Mo 


87r\/3 


are the renormalized coupling constant and mass, and Mq is the ultraviolet cutoh in units 
of m. The classical solution falls oh at inhnity as implying that a one-particle state 

has a mass rh. Thus, m is the mass renormalized on the mass-shell. 

Note that the second term in Eq.( p.39|) can be represented as 

- - —(140) 

2,2'Ky/^m? dz‘^ 

Thus, the change in the wall shape corresponds to a second order ehect in the shift of the 
wall center. 

The presence of the corrections to the shape of the domain wall has an impact on the low 
energy dynamics of the helds on the wall. The operator 0^0^ (z) is known to have a zero 


mode, corresponding to the breaking of the translational symmetry by the domain wall 


Exactly one bosonic and one (two component) fermionic zero modes exist, which rehects 
the balance of degrees of freedom in SUSY theories (the operator O^O^(z) does not have 
zero modes). This provides a zero mode superheld in the Lagrangian (|133|) denoted by $ 0 . 
This superheld (see Eqs. (P4|) , (p5|) ) does not have self-interaction, if all non-zero modes 
are omitted the Lagrangian for the zero mode superheld has the form 

1 


^ZM — j d^xcP‘d-{D^Qy 


(141) 


The zero mode, however, can mix with the non-zero modes through trilinear coupling and, 
thus, acquire a mass term induced by loop corrections with the massive particles inside 
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the loop. This is what is to be expected since we established that the shape of wall does 
experience corrections. The zero mode, being obtained by applying the broken symmetry 
generator {dz in our case) to the domain wall solution, also receives corrections. After we 
incorporate the loop effects into the effective Lagrangian, it is a certain linear combina¬ 
tion of all modes that remains massless and represents the Goldstone mode of the broken 
translational invariance. 


V. APPLICATIONS 

In this section we will consider some obvious applications of the ideas that were discussed 
above. 


A. Tunneling in non-supersymmetric non-relativistic quantum mechanics 

One of the most interesting results is the exact expression for e which can be obtained 
even in the absence of the explicit solution for the domain wall. This technique can be 
used in complicated multidimensional quantal problems with tunneling, for evaluating the 
tunneling exponents, in the situations where the WKB methods are inefficient, and the least 
action trajectories are not obvious. It is well-known, that in these cases one has usually 
resort to numerical solutions, which are quite complicated - since the effect sought for is 
exponentially small, achieving high accuracy is a non-trivial task. 

Exploiting the methods developed here the tunneling exponential in conventional (non- 
supersymmetric) quantum mechanics can be analytically obtained even without exact knowl¬ 
edge of the tunneling trajectory. 

Below we will consider a specihc example. It is quite clear however, that the situation 
is general. Consider four-dimensional (non-supersymmetric) quantum mechanics with the 
potential 

1/ = I - a + 7(X - $)2|2 + 1/9$ - - $)2|2 (142) 


45 



where the four degrees of freedom in question are the real and imaginary parts of X and 
All particle masses are put to 2. This is a reduction of the model discussed in Sect. III.B, 
see Eq. (pS]). The potential is horrible, it has two degenerate minima where the classical 
potential energy vanishes. We are interested in the probability of the system to tunnel 
from one minimum to another. As well-known, this probability is equal to the action of 
the classical trajectory connecting two minima in the Euclidean time. After the Euclidean 
rotation the equation of motion exactly coincides with that for the wall, what was in the 
case of the wall becomes the Euclidean time, and what was £ becomes the classical action. 
Since the solution is BPS-saturated the value of the action can be immediately read off from 
Eq. (1^ , without any calculations. The tunneling probability is 

oc exp(-| 8 a 3 / 23 -i^-i/ 2 |) _ 


B. Multiparticle production amplitudes at the kinematical threshold 


Non-vanishing corrections to the shape of the domain wall (Eq.( p.39|) ) have fenomeno- 
logical consequences. As was shown in Ref. m , the multiparticle production amplitudes at 
the kinematical threshold are related to the shape of the classical solution existing in the 
theory. We briefly review the main arguments leaving all the details aside. 

In the presence of the external source p{x) the multiparticle threshold amplitude at the 
tree-level can be written as 


An = {n\(j){x)\0) =f[f id‘^Xa)e ^ |p=o (143) 

where (jyd is the solution of the classical equations of motion in the presence of the external 
source p{t) and we made use of the reduction formula for the n-particle state with zero 
spatial momenta Pa = (a;, 0, 0, 0,). Since we are dealing with the threshold amplitudes, both 
p and (f) can be chosen to depend upon time t only, i.e. p = p{t) and 0 = 4’{t). The key point 
is that (pci can be shown to have the following expansion in powers of coupling constant A 
(cf. Eq.(|0l)) 
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( 144 ) 


4’ci = j(l) [l + X]c„(A z(t)Y 

\ n 


where 


z{t) = 


P 


(145) 


a;^ — Am? 

Factor 4m^ appears, since with our choice of the superpotential (|104|) , fluctuations over the 
vacuum have the mass 2m. Taking the limit as p —> 0 ,a; —> 2m and keeping z{t) hnite we 
obtain the threshold amplitude in the form 


(n|0(x)|O) = 


d 


dz{t) 


(i)d{z{t)) 1 ^= 0 , 


(146) 


and now is the solution of the (static) classical equations of motion without the source 


rf 

^Mt) = - (Wf iMt)), 


(147) 


with the boundary condition 


(j^ciit) —^ z{t) as A —> 0. 


(148) 


The solution of Eq. ( |147| ) can be obtained from the domain wall solution by substituting 
—A it. Rewriting it in terms of z{t) with the proper boundary conditions we obtain [|^] 

m 1 + Xz(t)/{2m) 


4’ci — 


(149) 


A 1 — Xz(t) / {2m) 

The loop corrections to the amplitudes An are, in turn, given by the quantum corrections 
to the shape of the classical solution 


, which are given in Eq . (11391) . All the details 


of the calculation can be found in Refs. |1l^ , p!3| , where non-supersymmetric theories were 
considered. Here we merely present the hnal result 


Ar, = (n|(/)(0)|0) = n\ 


X 


n—1 


1 + n{n — 1)- 


A2 A 


(150) 


2m) \ Stta/S/ 

Note, that the quantum corrections to the amplitudes An, are by factor 1/3 smaller than 
ones in non-supersymmetric theory []^. The result was known to M. Voloshin (private 
communication). 
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VI. “KINK” IN 1 PLUS 1 DIMENSIONS 


Now we turn to another example of supersymmetric theory, with SUSY being partially 
broken in the soliton sector. The two-dimensional analog of the Wess-Zumino model with 
the minimal {N = 1) supersymmetry was analyzed long ago The model is similar to 
the four-dimensional example. The theory possesses a classical solution - the “kink”, that 
breaks 1/2 of supersymmetry. However, contrary to the domain walls in four dimensions, 
the mass of the “kink” does receive a perturbative correction. 

The Lagrangian of the two-dimensional analog of the Wess-Zumino model describes the 
interaction of a real scalar held 0 and a real two component spinor T. It reads, in compo¬ 
nents, 

C = J - ^>V"(0)M| (151) 


where W(0) is the same as in Eq.( p.04|) and 


/ = 


a 


7 ^ = ia^. 


The consideration goes essentially along the same lines as in the previous section: we ex¬ 
pand the Lagrangian around the classical background and investigate the dynamics of the 
quantum huctuations. There are few changes, however. Now, there are only two super¬ 


symmetry generators Qi and Q 2 - The “kink” solution (|109|) breaks 1/2 of SUSY |^, and 
only the generator Q 2 survives. The surviving generator annihilates the “kink” solution, 
while applying the broken generator Qi to the “kink” will produce a fermionic zero mode, 
thus, transforming the “kink” into the left-moving fermion, and the “antikink” into the 
right-moving one. 

Let us note, that in contrast to the four-dimensional case, supersymmetry of the theory 
is not spontaneously broken. The “kink” solution has hnite energy and can be created from 
the vacuum and destroyed as a result of a held huctuation. In this sense the “kink” solution 
is a one-particle state, rather than the vacuum state of the theory, as it is the case with 
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the domain wall solution of the previous section. The “kink” can be pair-produced. Super- 
symmetry algebra will be realized linearly on the states with different number of “kinks” 
and “antikinks”. We are interested, however, not in the general supersymmetry of the the¬ 
ory, but in the question how it is realized in the given soliton background. In the given 
background Qi must be considered as spontaneously broken. 

The question of the corrections to the mass of the “kink” was addresses more than once 
in the literature (see |^, |^, and references therein). There are two basic approaches that 
differ in the way the boundary conditions on the quantum fluctuations around the “kink” 
solution are imposed. 

In the hrst approach, that was taken in [^, non-supersymmetric boundary conditions 
were imposed at spatial inhnities. Since supersymmetry is explicitly broken by the choice 
of the boundary condition, the volume effects associated with the production of virtual 
fermions and bosons will introduce a non-vanishing mass correction. 

To see this let us look at the following eigenvalue equations (cf. Eq. (|121| ), ( p.22|) ) that 
can be easily obtained by expanding (|151|) around the “kink” up to the quadratic terms : 


0-‘0\z)n = AfS-i , O^OHz)^- = (162) 

where 

(Note that are real one-component helds in contrast to the helds dehned in the pre¬ 
vious section that had one independent complex component. With our choice of 7 matrices 
= Ti and T+ = ^2 ■) 
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FIG. 2. Eigenvalues of two operators. 

If non-supersymmetric boundary conditions 

0(-L) = (^(L) • (153) 

<P'{-L)=<P'{L) 

are imposed (L —> cxo), then in the limit of inhnite volume two operators 0^0^ and 0^0^ 
will have different density of the eigenstates. This difference can be readily deduced from 
the asymptotic form of the eigenfunctions of two operators and is equal to 


d n~^ dn 


(154) 


d k d k 7r4 + fc^ 
where k is the spatial momentum in units of m. Integrating over the spatial momentum we 
obtain, 


N+ -N- = I 


(155) 


where {N-) is the number of the continuous eigenstates of the operator 0^0^ (or 0^0“^). 
This difference in the density of the fermionic and bosonic eigenstates will produce a non-zero 


mass correction that was derived in 31 


An alternative logical step would be to impose a supersymmetric boundary conditions. 
Then, the consideration will run in parallel to what was done with the domain walls in 
the previous sections. This approach was undertaken in where the supersymmetric 
boundary conditions suggested in were modihed to be 
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( 156 ) 


0(-L) = 0(L) 

same as before for fermions 

Then, all volume contributions in the one-loop correction cancel out, the boson correction is 
canceled by the fermion one. Mass renormalization, nevertheless, still occurs in [Q due to 
a surface term that appears after integrating the term in Lagrangian (|151|) by parts 


and using (|156| ). 

Let us note, that the cancelation of the bulk effects is usual for the theories with un¬ 
broken supersymmetry, the same situation occurred when domain walls were considered. 
However, the appearance of the boundary correction is specihc for supersymmetric theories 
in two dimensions, where superspace can not be divided into chiral subspaces. In four di¬ 
mensions, quantum corrections can not induce non-vanishing (<!>"') for the chiral held $ at 
the boundary, hence no boundary correction can appear (in other words the superpotential 
is not renormalized in four dimensions). In contrary to four-dimensional case, in two dimen¬ 
sional models with the minimal supersymmetry no chiral held exists, and, thus, (<h”) can 
have non-zero value, inducing non-vanishing boundary terms. In two-dimensional models 
with the extended supersymmetry [N = 2) the quantum corrections from the boundaries are 
absent much in the same way as in the domain wall problem, and the soliton mass remains 
unrenormalized. 

There is a subtlety in imposing the supersymmetric boundary conditions. The point is 
that the soliton mass correction is given by the surface term, i.e. it is due to the region most 
ahected by the choice of the boundary condition. This is true for any supersymmetric choice 
of the boundary conditions. The answer is saturated on the boundaries and, thus, largely 
affected by particular choices of the boundary conditions. It is necessary that the given set 
of boundary conditions corresponds to a pure soliton state rather than some excited soliton 
state. The following simple argument illustrates that this is not the case with the boundary 
conditions (|156|). Indeed, Eq. (|156|) implies that the expectation value of the operator (TT) 
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is the same at ±cx). On the other hand, since the fermion mass changes its sign as one goes 
from plus to minus inhnity, the expectation value of in the pure soliton sector should 
also change sign. This inconsistency implies, that the boundary conditions 


project out 

an excited soliton state containing a fermion pair. Physically, it is clear that the problem 
of calculation of the soliton mass must be formulated in a way independent of the specihc 
choice of the boundary conditions. This can be done by considering a local form of the BPS 


saturation condition. For more details see 39 


Let us note, that in line with the arguments presented above, the subtlety with the 
boundary conditions is immaterial in the four- dimensional theories, and, thus, need not be 
considered when dealing with the domain walls. In four dimensions the expectation value of 
the operators (<!>"') is zero due to the same argument that leads to the non-renormalization 
of the superpotential in four dimensions. Thus, the contamination of the soliton does not 
occur, and the particular choice of the boundary conditions does not affect the solution. 


VII. CONCLUSIONS 

The existence of the central extension in the N = 1 superalgebra in four-dimensional 
gauge theories was hrst observed in Ref. [0,|1. We present a detailed derivation of this exten¬ 
sion in the general case: non-Abelian gauge theories with a superpotential. Our derivation 
takes advantage of the well-established “geometric” anomalies, such as the anomaly in the 
trace of the energy-momentum tensor. In fact, we show that the latter is in one-to-one 
correspondence with the central extension in the N = 1 superalgebra which we thoroughly 
discuss here. 

In few cases when the diagrammatic computation is simple, the result agrees with the 
general answer obtained indirectly. The direct diagrammatic derivation for the general case 
is the subject of further work. 

The fact that the superalgebra is extended and the central charge is non-zero leads to 
far-reaching consequences: the existence of the BPS-saturated topological defects. In-depth 
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studies of the BPS-saturated domain walls is carried out. We derive the creek equations and 
hnd necessary conditions for the solutions to exist, and other useful criteria. BPS saturation 
means that solutions partially preserve supersymmetry, 1/2 of SUSY remains unbroken in 
the given background. This, in turn, entails exact calculability of the domain wall energy 
density. It is related to the central charge for the corresponding topological defect. The 
domain wall energy density receives no quantum corrections. 

The shape of the domain walls can not be exactly deduced from the central extension and 
is subject to perturbative loop corrections even in the presence of residual supersymmetry. 

We discuss in great detail the generalized Wess-Zumino models. Apart from the domain 
walls, string-like topological defects, that are BPS-saturated, are shown to exist. This is 
a novel step, since the BPS-saturated strings were not discussed in Ref. [0]. It is clear 
that consideration of the BPS-saturated topological solitons in the supersymmetric {N = 1) 
Yang-Mills theories only begins, and will be continued further, see also Ref. . 


In the minimal Wess-Zumino model, where the solutions of the classical equations of 
motion are known analytically, all conclusions based on the analysis of the central extension 
can be explicitly checked. In the strongly interacting gauge theories, to explicitly reveal 
the way the central extension works, one has to resort to effective Lagrangians of the type 


suggested by Veneziano and Yankielowicz, as was done in |^. We do not rule out that 
other, more elegant methods will eventually emerge. 

Note Added. A few days ago, when the write up of this paper was being completed, 
Witten’s work was submitted to hep-th, where constructions similar to those discussed 
above were discussed in the context of the D-brane-based approach to QCD. 
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APPENDIX A: WESS-ZUMINO MODEL AND SUPERSYMMETRIC 
GLUODYNAMICS IN SPINORIAL FORMALISM 

1. Notation 

We hrst specify our notations and conventions. They are close but not identical to these 

of Bagger and Wess |Q. The main distinction is that we choose the metric to be (H-). 

There are also distinctions in normalization, see Eq. ( |A190 . The left-handed spinor is denoted 
by undotted indices, e.g. r]p. The right-handed spinor is denoted by dotted indices, e.g. 
(This convention is standard in supersymmetry but is opposite to one accepted in the text¬ 
book 1^). The Dirac spinor 4/ then takes the form 


T = 




(Al) 


Lowering and raising of the spinorial indices is done by multiplying by the Levi-Givita tensor 
from the left. 


X" = e“^X/3 , Xa = 


(A2) 


and the same for the dotted indices, where 


_ _g/3« ^ _ _ eu = 1. 


(A3) 


The products of the undotted and dotted spinors are dehned as follows: 


hX = h"Xa = -haX" , VX = haX“- 


(A4) 
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Under this convention {r]x)^ = XV- Moreover 

IP0f> = . 

(A5) 

The vector quantities (representation , i)) are obtained in the spinorial formulation by 
multiplication by 

( )a/3 = { 1 > ^ }a/3 (^6) 

where r stands for the Pauli matrices, for instance, 

(A7) 

Note that 

A.A'-. (A8) 

The square of the four-vector is understood as 

= \a^pA^^. (A9) 

If the matrix is “right-handed” it is convenient to introduce its “left-handed” coun¬ 

terpart, 

(B- f- = { 1 , -f (AlO) 

The matrices that appear in dealing with representations (1 , 0) and (0 , 1) are 

{^Tp = Ta/3 , (All) 

and the same for the dotted indices. The matrices (cf)"^ are symmetric, (cf)"^ = ((t)^“. In 
explicit form 

[ayP = { , -fl , -r^ }ap-, , fl , -r^ }ap . (A12) 
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Note that with our dehnition 


}«/3; (<?)«/? = { -'T^ , . 

The left (right) coordinates xl^r and covariant derivatives are 


^aa 2 i 6^6^ i (2^/?)ao Xaa T 2 i da^a i 

^ ~ * ^ + * 6^ (9q,q,, 


SO that 


{Do, D^} = 2ido,a , 

Dpi^Xl{)(xa 0 , Djj{^XR^(xa 0 , 

Dp{,^Fi)eta 4 i ^a^0a > Dp{^Xi{)ota 4 i (^a^pa- 

The law of the supertranslation is 

6 —> 6 -\- e , 6 —> 6 -\- e , 

^a0 ^ 2:^/3 - 2 i 9o,e0 + 2 i ej^ , 

(3;i)a/3 ^ i^L)ot0 - 4 ^ , 

i.^R)ap (^R)a0 + 4 * 

It corresponds to the inhnitesimal transformation of the superfield in the form 

5$ = i(eQ + eg) $ 

where 

Qa = 4 " ^atj^^i 

The integrals over the Grassmann variable are normalized as follows 

J (Pe6^ = 2, J d^e= A, 

and we dehne 

{...u = i /rfM-..}. {...}d = \ jcf0{...}. 


(A13) 

(AM) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 
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2. Wess-Zumino model 


The minimal Wess-Zumino model describes (self) interaction of the chiral superfield 


^{xl, 9) = (j){xL) + V29'tp{xL) + 9‘^F. 


(A20) 


The action is given in Eq. (P). In components it has the form, 


C = + f+f 

+ {f w\^) - + H.c\. 


(A21) 


The model possesses a supermultiplet of currents, 


Jaa = ^ {T>a<I>T>a<h - 


(A22) 


The 9 component of the supermultiplet is related (cf. Eq. (A38|) ) to the supercurrent, that 
has the form 


4/3/3 = 2^2 { [(440+)3/’/3 - i e/3a^4 

1 r 


6 1 


4/3(V’/30 ) + 4/3(V'a0^) - 


(A23) 


where 


F = 


dW 

90+ ’ 


and W is the superpotential. In the spin-vector form the supercurrent is 


(A24) 




a/3/3- 


(A25) 


This current is conserved. 


d, K ^ = 0, 


(A26) 


The supercharge is defined 


Qa = J d^xJl- 


(A27) 
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Terms in the second sqnare brackets are the total spatial derivatives and, thus, do not 
contribute to the supercharge. These terms have to be added if we require the supercurrent 
to enter the same supermultiplet as the energy-momentum tensor. Then at the tree level, the 
current has vanishing trace in the the theory, that is classically conformally invariant. 
The relative coefficient in front of the second bracket is most readily hxed by the requirement 
that reduces to the equation of motion. 

The anticommutator of Qa and is standard, 

{QaQp} = 2P^^, (A28) 

while that of Qa and Qp contains a central extension 

{Qc Q^} = I |w - . (A29) 

At the loop level, due to the quantum anomaly, an additional term on the right-hand side 
appears; being combined with (l/3)d>9W/95), this additional term is a total super derivative, 
that gives no contribution to the central charge. Hence, we arrive at the following exact 
relation 

{Qa Qp} = (-4*)(a)a/5 J d^xVW. (A30) 


3. Supersymmetric gluodynamics 


This model (as well as its A^ = 2 super generalization) was hrst introduced by Ferrara 
and Zumino 0. The superheld Wa which includes the gluon strength tensor is 

1 


Wa = 


D"’- (e-'' e'') 


(A31) 


where V is the vector superheld; in Wess-Zumino gauge 


i/ = -29^9^Aaa - 2i9\9X) + 2i9\9\) + 9^9^D, 


(A32) 
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V = and T“ stands for the generators of the gauge group G, 

Tr ^ 5“^ 

for the fundamental representation. If the matter held is denoted by S, the supergauge 
transformation has the form 

S , ^ ^ (A33) 

where A is an arbitrary chiral superheld (A is antichiral). In components 

= i(\a + iOaD - (A34) 

where Aq, is the gluino (Weyl) held, Vaa is the covariant derivative, and Gap is the gluon 
held strength tensor in the spinorial notation. 

The standard gluon held strength tensor transforms as (1, 0) + (0,1) with respect to the 
Lorentz group. Projecting out pure (1, 0) is achieved by virtue of the matrices, 

= (J - iB){SU. (A35) 

Then 

G^f^Gap = 2{B^-E^+ 2 i EB) = G^M^^ - i 

where 

G^u = -^^uapG'^^ , (eoi23 = ~1)- (A36) 

The Lagrangian of the supersymmetric gluodynamics is given by 

where d is the vacuum angle. 

The supercurrent superheld takes the form (the trace is taken over the color indices only) 

Jad = -^ Tr (iT^e-'" Wa e^) = 

-^A“A“ - i (^Jpaa - lepae^^Js,a) + H.G.^ - 

0^0^ [Ja6.pp - +■■■ (A38) 
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where Jpaa is the supercurrent 


4^ c — 'j 

Jl3aa = (A39) 

and energy-momentum tensor, 

'^aaPP ~ ~ * pp^a “ i {T^ (A40) 

The energy-momentum tensor in the spinorial notation is related to 9^y in the following way 

{»" + »“9” - (A41) 


where i, j, fc = 1, 2, 3 , g^i, is the metric tensor and matrices {(T'^)ap are dehned in Eq. (|All|) . 

The supercurrents above are normalized in such a way that the standard commutation 
relation holds, 


{Qa Qa} 2 Petal 


and the central extension of the superalgebra then is 

-iN 


{Q» Qi,} = (yx-) 


(A42) 


modulo full superderivatives. Finally, we give (for future applications) a few expressions for 
SQCD (i.e. supersymmetric gluodynamics plus matter). 

Consider for simplicity SU{2) gauge group and one matter held in the fundamental 
representation (two subhavors). The corresponding chiral held carries the color and subhavor 
indices. The Lagrangian is 
1 


Cm = 


-J J / cCO S^^Saf + H.c)j 


(A43) 


In components it has the form 


C = 




+ 




iV2{(l)'l\'lpi 0jA-?/>2 + 02AV’2) + 

(0+T>i-F 0+T“02) + 

m ((^1^2 -F 0+^2+ -h 02^1 -F + '0lV’2) 


(A44) 
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The supercurrent is now modified, 


>^0/3/3 “ ^ ~ 3e/3„T)“A^) + 


^ f(^a/ 3 ‘^^)V ’/3 - i ^PaF'ljji 


and the central extension is 


v^r 


^a/3(V'/30^) + ^/3/3(V’a0^) “ 


(A45) 


{QaQ0] = (-4*)(C?)„^ / d^xV 


W - f^ TrlT^ 

lOTT^ 


J 3 e=0,6»=0 


(A46) 


It is also useful to write down the Konishi anomaly equation in the superfield formalism: 


^2 gaf\^ ^ 4 ^ gaf ^ ^ 


(A47) 


For future references, let us also give the expression for the Lagrangian of the SQCD 
with the arbitrary gauge group and superpotential W 

/:=|^Tr Jd^e +B.c.j + ^Jd^0 ^Q.eyQ, + Qjd^0 W(Qi)+ff.c.j. (A48) 


APPENDIX B: 


To obtain the three-dimensional Lagrangian we substitute expansion (|123|) into La¬ 
grangian (|120|, |136|) and integrate over 2;. We will need to compute the following set of 
integrals (brackets stay for the integration over z), 


( '} •} { '} •} ( ^k )? 

( bj Ojj 1 ( 6; ) . . . 


(Bl) 


where a* and bi are the eigenmodes (|124|) and (pu, is the classical solution. Fortunately, to 
establish the existence of SUSY we do not have to compute the integrals, but rather to 
establish certain relations between them. With the aid of the completeness relation 


\n an{z) ttniy) = S{z - y) 


(B2) 
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where A„ is the normalization factor, and repeated use of Eq.( p.24|) all the integrals can be 
reduced to the following two: 


^ijk \J Afc ( Qjj -^ijk ^k ( ^k )• 

Coefficient Aijk is symmetric, while Dijk is symmetric with respect to the last two indices. 
Note also, that since the operator 0^0^(z) does not have a zero mode, then 


-Dooo — ^000 — 0. 


(B4) 


With the aid of the following redehnitions: 


— / t : ? Hn — 


V2 


dn bji 

V2 ’ 

Bijk Dijf^ Djij^ Dj^ij 


~ H“ byi ~ 

1 — 


Cijk Dij]^ Djij^ Dj^ij A 


ijk- 


(B5) 


the Lagrangian (|120| , |136| ) can be seen as a component expansion of supersymmetric La- 
grangian ( p.33|) . 


APPENDIX C: PROPERTIES OF THE OPERATORS O'^O^ AND O^O'^ 

The operators 0‘^0^{z), 0^0^(z) have the following form: 

O^O^(z) = -dl + >V"2(0^) - WW'\K). 

0‘^0\z) = -dl + + WW"(0^). (Cl) 

After we substitute the classical solution 

Tfl 

(pw = —tanh(m^), (C2) 

A 

we arrive at 
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= -dl + 4m2 


2w? 


0^0\z) = -d^ + 4m^ 


cosh^(m2;) ’ 
Qm? 

cosh‘s (mz) 


(C3) 


Operator 0‘^0^{z) has two discrete eigenvalues and continuum, 


Aq =0; A]^ = a/Stu; A^ = a/Zc^ + 4m^. 


(C4) 


Operator 0^0^{z) has only one discrete mode and continuum, 


A;^ = VSm; Xj^ = \/kX + 4m?. 


(05) 


It is obvious that if bn is an eigenstate of 0^0‘^{z) then an = 0‘^bn is an eigenstate of 
0^0^ {z) with the same eigenvalue. Thus, two operators have exactly the same spectra, 
except for the zero mode of 0‘^0^{z) . This is the type of symmetry one encounters in 
supersymmetric Quantum Mechanics p8 . 
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